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We develop a linear response formalism for nuclear spin diffusion in a dipolar coupled solid. The 
theory applies to the high-temperature, long-wavelength regime studied in the recent experiments of 
Boutis et al. [Phys. Rev. Lett. 92, 137201 (2004)], which provided direct measurement of interspin 
energy diffusion in such a system. A systematic expansion of Kubo's formula in the flip-flop term 
of the Hamiltonian is used to calculate the diffusion coefficients. We show that this approach is 
equivalent to the method of Lowe and Gade [Phys. Rev. 156, 817 (1967)] and Kaplan [Phys. Rev. 
B 2, 4578 (1970)], but has several calculational and conceptual advantages. Although the lowest 
orders in this expansion agree with the experimental results for magnetization diffusion, this is not 
the case for energy diffusion. Possible reasons for this disparity are suggested. 
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I. INTRODUCTION 

Spin diffusion occurs for conserved quantities in para- 
magnetic spin systems. For dipolar-coupled nuclear 
spins in high external magnetic fields, these quantities 
are the interaction energy and the spin component along 
the field. In a solid, diffusive behavior arises from energy- 
conserving two-spin flips of anti-aligned spins, which pro- 
duce a dynamics analogous to a random walk.-'^ Recent 
experiments have provided the first direct measurements 
of magnetization^ and energy'* diffusion in a dipolar- 
coupled spin system. These experiments were carried 
out on very pure single crystal samples of the insula- 
tor calcium fluoride, a substance that has been studied 
extensively by the NMR community for the past sixty 
years, and whose structural properties are well knownj^ 
Calcium fluoride is an ideal material to study because 
of its very long spin-lattice relaxation time (Ti ^ 100 
s at room temperature in these experiments), its simple 
structure (the fluorine sublattice is simple cubic), and 
the high accuracy to which the spin-1/2 nuclei behave as 
ideal magnetic dipoles. Since the timescale on which dif- 
fusion was observed (on the order of seconds) was much 
shorter than Ti , the spin system was well isolated and the 
evolution coherent. This is the regime in which Waugh 
and co-workers demonstrated that spin diffusion is re- 
versible (by a sign change of an effective Hamiltonian via 
multiple-pulse NMR methods) Nevertheless, the trans- 
port of magnetization in three dimensions appears indis- 
tinguishable from an incoherent process, as many calcu- 
lations and experiments have shown i2i2ii24iiii24iiiiiiSii& 
In contrast, we find that this does not seem to hold for 
energy transport, where despite diffusive behavior the use 
of thermally averaged correlation functions leads to re- 
sults that are inconsistent with the experiments of Boutis 
et al. 

This paper is organized as follows. In Sees. II and 
III we outline the linear response formalism for spin 
diffusion in a dipolar-coupled solid, including a deriva- 
tion of the energy and magnetization current operators. 



This formalism applies to the high-temperature, long- 
wavelength regime studied experimentally. It is equiva- 
lent to the density matrix approach of Lowe and Gade,— 
and Kaplani4 (LGK), as we prove in the Appendix. Our 
formulation has the advantages that the long-wavelength 
limit is built-in and the application to inter-spin energy 
is straightforward. In Sec. IV, we derive the expansion of 
the Kubo formula in powers of the flip-flop term of the 
dipolar Hamiltonian, which we use for numerical evalua- 
tion of the diffusion coefficients. The diffusion coefficients 
of magnetization and inter-spin energy obtained from this 
expansion to two leading orders in the flip-flop are given 
in Sec. V. Their numerical values are calculated in Sec. 
VI, along with an estimate of the errors. Our main find- 
ing is that, although the expansion for the magnetization 
diffusion coefficient reproduces the experimental results, 
the series for inter-spin energy does not. In Sec. VII 
we discuss this result and comment on how our method 
complements previous approaches to the problem. 



II. KUBO FORMULA 

The linear response diffusion coefficients of a continu- 
ous system may be obtained from Kubo's formulajii 

/p°°dt/d3r/dV(j,(r,t)j,(r',0)),, 

/d3r/d3r'(5(r,0)5(r',0))e, ' ^ ^ 

where iS'(r) is the operator representing the local den- 
sity of a globally conserved quantity. In our case this 
is either the energy {S = Ti) or the component of spin 
magnetization along an external field {S = M), and jz(r) 
is the corresponding current density along the transport 
direction. The angular brackets denote averaging over an 
equilibrium statistical ensemble at the appropriate tem- 
perature, which is infinite for the case we consider. The 
discrete version of this formula is more directly applica- 
ble to our problem since the spins reside on a lattice, and 
is given by replacing r with a lattice site index. However, 
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the continuum representation is more useful for deriva- 
tions, and is completely equivalent. 

Eq. ^ is the correct formula for the diffusion coeffi- 
cient of the quantity S if the following assumptions hold. 

(1) The correlation function of 5* is known to have a 

diffusive form (i.e. a diffusive pole). 

(2) The system may at all times be described by a sta- 

tistical ensemble that is sufficiently close to equi- 
librium (linear response regime). 

Assumption (1) may be checked experimentally, and has 
been verified for both spin-spin and energy-energy corre- 
lators of dipolar-coupled spins in a solid. '^''^ Assumption 
(2) is more subtle, as it rests on the validity of the ergodic 
hypothesis, which has received much attention recently 
in the context of lattice spin systems in dimensions 1 - 
3j18j19j20j21j2^ As discussed below, we believe that it may 
not be valid for the diffusion of energy. 

Experimentally, spin diffusion has been measured by 
observing the relaxation of initial states varying sinu- 
soidally in real space with a given wavevector Such 
spatially-inhomogeneous states are represented mathe- 
matically as perturbations on the infinite temperature 
equilibrium state, poo — 1/2^, where N is the number 
of spins. We have 



p(0) = Poo + 6piO). 



Two possibilities exist for Sp(fi), corresponding to long- 
wavelength fluctuations in the two conserved quantities, 

SpM ik,0) = e J d^rcos(k • r)M(r), (3) 

SpniKO) = e J rf^rcos(k- r)H(r), (4) 

where e is a small quantity of order jhBo/2^kBT. Here 
Bo is the external field, and 7 is the gyromagnetic ratio of 
the nuclear species of interest (7 ~ 2.51 x 10'* rad Hz/Oe 
for *^F). Using typical experimental values of Bq — 
1 T, T = 300 K, we estimate e ~ 10'^ - 10"^. The 
important length scales are related by L >> fc~* >> a, 
where L is the sample size and a is the lattice spac- 
ing, a = 2.73 X 10~®cm for the Fluorines in CaF2, and 
L ~ 0.1cm, ^ 10~^cm in the experiments. The re- 
alization of such states is discussed in detail by Boutis 
et alX In practice, k is parallel to the external magnetic 
field, and transport is measured in the same direction. 
As shown in the Appendix, expectation values taken with 
respect to the states in Eqs. and I0J give rise to the 
equilibrium-averaged correlation functions appearing in 
Eq. iPl. 

We find it plausible that assumption (2) is valid when 
the spatial profile of magnetization or energy in the ini- 
tial state varies sufficiently slowly, and when this state 
has no long range correlations. While this is possible for 
magnetization, the creation of spatial inhomogeneities in 



the interaction energy requires NMR techniques that in- 
troduce short - ranged correlations between the spins. 
Some of these correlations persist due to energy conserva- 
tion, and these are the ones relevant to energy transport. 
Other correlations are assumed to decay rapidly and to 
be unobservable. A careful analysis of the initial state 
involved in the energy diffusion measurements is needed 
to determine whether assumption (2) is satisfied. On the 
other hand, the extent to which our results agree with 
experiment may be viewed as a constraint on its validity. 



III. CURRENT OPERATORS 

In order to use Eq. ^ we need to obtain expressions for 
the current operators jz(r). We start with the rotating 
- frame Hamiltonian of interacting magnetic dipoles in a 
large external magnetic field iSiS^ 



1-1= ^ Hjj, 



(5) 



'H^J = hj - - [h+Ij- + Ij + h-)j ■ (6) 

The latin indices run over all lattice sites and the 
are operators for spin-1/2, satisfying the commutation 
relations [liajljjs] = Sijli-y, where a, 7 is any cyclic 
(2) permutation of x, y, z. The Ij± = Ijx ± iljy are raising 
and lowering operators. The coupling constants are given 

by 



-f'^Ti^ 1 - 3fi 



(7) 



where is the displacement vector between lattice sites 
i and j, and unit vectors are denoted by a caret. We take 
the z-direction to be along the external magnetic field. 

The total z-component of spin and total energy are 
conserved quantities of the Hamiltonian, Eqs. © and 
©. This implies that they satisfy a local continuity 
equation^ We define spin and interaction energy den- 
sities that depend on continuous variables fi2i24 



/c(r,t) 



6{r ~r.,)ha{t), a = z,+,- 



(8) 



(9) 



The operators in these equations and throughout this 
section are in the Heisenberg representation. 

The densities of magnetization and interaction energy 
are given in terms of the spin densities as 



M{r) = jhh,{r), 

n{r)= j d'v'h{v-v'){h{v)h{v') 



(10) 



[/+(r)/4r')+/+(r')/-(r)]}, (11) 
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where b{r — r') = ^ ■ j 6{r — ri)S{r' — Yj)hij. We avoid 
writing the time dependence of the operators where there 
is no possibihty of confusion. Next, we use the Heisen- 
berg equation of motion, 



dM{r,t) _ i 

dt ~ ~h 

dn{r,t) _ i 

dt ~ ~h 

to obtain 

dM{r) _ i^h 



Mir,t),n], 



(12) 
(13) 



dt 



^ J d3r'6(r-r')(/+(r)/_(r') 

-/+(r')/-(r)) 



(14) 



dnjr) 
dt 



= (fiv'(Pr"b{Y' - r)b{r" - r') 

x{4(r')(/+(r)/_(r")-/+(r")/-(r)) 
+2/.(r) (/+(r')/-(r")-/+(r")/-(r')) 
+2/.(r") (/+(r)/_(r') - /+(r')/- (r))}. (15) 

In order to write Eqs. (|14|l and H15|l as continuity equa- 
tions in the usual form, they may be integrated over an 
arbitrary volume and the result expressed in terms of a 
surface integral, whose integrand is the current density. 
This gives the current density operators, 

j(M)(i.,t) - ^ y dh'bir - r')(r - r') 

x[/+(r)/_(r')-/+(r')/-(r)], (16) 

j(^)(r,t) - ^ j d'v'SY"h{v' - v)h{v" - r')(r - r") 

x{/4r')[/+(r)/-(r")-/+(r")/-(r)] 
+2/,(r) [/+(r')/-(r")-/+(r")/-(r')] 
+2W) [/+(r)/_(r') - /+(r')/-(r)]}, (17) 

where the superscripts M and TL denote magnetization 
and energy current, respectively. The continuity equa- 
tions then take their usual form, 

^^+V.j(^)(r,i)=0, S = M^n. (18) 

Similar results have been obtained by Furman and Goren 
by a different method 

IV. PERTURBATION THEORY 

Using an interaction representation introduced by 
Lowe and Norberg,^^ we can expand Eq. in powers of 
the flip-flop term of the Hamiltonian. Following Ref. 13 
we define 



(19) 



— bij lizljz, 

n2^-\ b,j{h+h-+U-Ij+). (20) 



Using the notation A{t) = g-'^i*/'* A(0)e''^i*/'"\ we may 
write any operator in the Heisenberg representation as 
an infinite series in 7^2, 



A{t) = e^-'l'^ \ ^(0) -y^ dt^ [7^2(ti),^(0)] 
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ti 



dti / dt2[n2{ti),[n2{t2),A{o)]] 



Jo 

iHit/n 



(21) 



Using this expansion for jz{t) in Eq. JQ) gives a pertur- 
bation series for the diffusion coefficient. 



D 



dt{jz{0)l{t)) 



dt dt,{[n2iti),jzmjzit)) 

n Jo Jo 

dt / dt, / dt2{[n2{ti),[n2{t2),jzm] 



(22) 



"'0 "'0 



<]z{t)) 



where the spatial integrations have been suppressed, and 
angular brackets denote an equilibrium average. The 
operators 7^2(0 ^^'^ Jzit) can be evaluated using the 
idcntity,iiS& 



L„(t), (23) 



where 



(24) 



It is advantageous to approximate the operator Lij (t) by 
a c-number, equal to its thermal average, Lij (t) ~ Gij (t). 
At infinite temperature, 

Gzj{t)^ Y[ cos{b,m - bjm)t. (25) 

The neglected q-number terms are expected to be ap- 
proximately 20 — 25% smaller, as discussed at the end of 
Sec. IV. 

It turns out that Gij (t) depends very weakly on its 
indices. In order to make the calculation of higher order 
terms in the perturbation series tractable, we therefore 
replace Gij[t) by G{t). G is a suitably defined average 
over all the G^-'s, as discussed in Sec. VI. The overall 
error introduced by approximating the higher order terms 
in this way is therefore small (we shall see in Sec. VI that 
it is less than 10%). The resulting perturbation series 
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contains only even-order terms, 



B. Energy Diffusion 



POO 

F= Git)dt, (27) 
Jo 



where [A,B]n = [yl, B]...]] is a commutator 

with A taken n times. It is worth noting that this se- 
ries may be written in a simple closed form, 



D 



(sW) 



(28) 



which shows that our system is approximately equivalent 
to one with purely flip-flop interaction, evolving for a 
finite time, F. 



V. ANALYTICAL RESULTS 

In this section we give analytic expressions for the dif- 
fusion coefficients to the two leading orders using Eq. H22|l 
and (|26|l . We are interested in the infinite temperature 
case. Therefore, (• • •)eg = tr{---}/2^, where N is the 
number of lattice sites. Each term in Eq. H26() then con- 
tains a trace over a product of spin operators at different 
lattice sites. The evaluation of such traces is a tedious 
but straightforward task. We do not discuss it further 
here. 



A. Magnetization Diffusion 

The denominator of Eq. ||T} for S* M at T = CX3 IS 
/ d^r J d^r' tr {M(r, 0) A/(r', 0)}/2^ = N-f^h^A. 

The lowest order term is calculated exactly, by insert- 
ing the expressions for the Hamiltonian, Eq. H19|l . and 
the magnetization current, Ea. (|16|l . into the first line of 
Eq. and applying the identity, Eq. (j2Sl)- We obtain 



D 



(0) 
M 



4 ^ 



(29) 
(30) 



where Fy is given by Eq. lf77|) with Gij{t) replacing G{t). 
The lattice site labelled is arbitrary, due to the trans- 
lational invariance of an infinite lattice. This result is 
identical to LGK's, as expected from the equivalence of 
the two methods (see Appendix). 

The next term, obtained from Eq. H26|l using all of the 
approximations discussed in the last section, is 



For S — Ti, the denominator is 
J d^rj d^r'tT{n{r,0)nir',0)} = (3/16) x 2^ x 
^"ij- The lowest order term is 



D. 



H 



48 ^i,) \i,j,fc(i#j#fe) 



—2^ BijkBiikRikji 



(32) 



where 



B'ijk — bijbjkZik ~t- 2bijbikZjk -\- 'ZbikbjkZij . 



Rkl-ij = 



dtKki-ijit), 



(33) 



Kki-ij{t) = -4:tr {IizIjzLki{t)}\i=£j^k=£i 



sin(AA;;;ii) sm{Aki-jt) J| cos(Afe;;pi),(34) 

p{P7^'i,j.k,l) 

Afci;p = bpk - bpi. (35) 

It is useful to have approximate analytic expressions 
for the integrals in Eq. H27|) and Eq. I|33|) . It has been 
shown that the saddle point approximation to Fij is quite 
accurateiiSiii It is equivalent to replacing in the in- 
tegrand by a Gaussian, 



Gy(i) « exp -i ^ {bk^ - bkjff 

V k(k^i,j) ) 



(36) 



To evaluate Rki-ij, we replace the product of cosines in 
Kki:ij (t) by a Gaussian, as above. Since this Gaussian 
cuts off the integral at times t ^ h/b, where b is of the 
order of the nearest neighbor coupling strength, we can 
expand the sine terms around t = 0, which yields 

Kki-ij{t) w {bik - bu)ibjk - b.ji)t^ 



X exp 



\ P(p¥'i,3,k,l) J 



These approximations give 



Rkl-ij -1 / — 



'^T,k{k^i,3)(bkt - bkjY 

{bik - bii){bjk - bji) 



2 / \ 3/2 ■ 

J2p(p^iJ,k,l}i^pk ~ Kl)^ 



(38) 
(39) 



D 



(2) 
M 



I 2 E 

i{i^O) 



2 1,2 2 



(31) 



Equation H38|l has been derived previously»iii4 Equa- 
tions H38|) and (|39|l allow a more rapid numerical evalua- 
tion of the expressions for the diffusion coefficients than 
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by numerical integration of Eqs. H25() and H34|l . By nu- 
merical integration on cubic lattices of between 5'^ and 
ll'^ spins, we have verified that Eq. (|38|l approximates 
the exact value with an error that is less than 2% and 
Eq. H39() approximates the exact value to within 10%. 

We found numerically that the sum over R terms in Eq. 
(|32|l was about 20 — 25% in magnitude of the sum over 
F terms, for both the (001) and (111) directions. The 
R terms arise from the lowest order (in J^^ operators) 
q- number correction to the approximation, Eq. H25|l . to 
Lij{t). We therefore expect this same number to also be 
a good estimate of the error in this approximation. 

Using Eq. H26|l for the next order term, we obtain 



TABLE L Summary of the theoretical and experimental re- 
sults for the spin diffusion rate of energy, Dh , and magnetiza- 
tion. Dm for a single crystal of calcium fluoride. Theoretical 
values have been obtained by numerically evaluating Eqs. 
and igOll, and Eqs. ^ and igB, 

using Eq. (1381 for Fij and 
Eq. 13911 for Rikmj , and using finite size scaling to extrapolate 
to the infinite lattice limit. 



^(2) _ 



F3 



(3!)(192)E 



/)2 



Bqul^^^uqbqlBqulBulq) 



uq q' 



u,q,l{u^qy^l) 



Theory 


[001] 


[1111 

L J 


Dqoi /Dill 




(x 10-12 cmVs) 


8.4 ±0.2 


7.9 ±0.2 




J)(2) 


(x 10-12 cm^/s) 


-0.3 ±0.1 


-0.5 ±0.2 




n(o+2) 


(xl0-i2cm2/s) 


8.1 ±0.3 


7.4 ±0.4 


1.1 ±0.1 




(xl0"i2cm2/s) 


19.8 ±0.7 


12.5 ±0.3 






(Xl0-i2cm2/s) 


-0.7 ±0.3 


-1.1 ±0.5 




n(o+2) 


(Xl0-i2cm2/s) 


19.1 ± 1.0 


11.4 ±0.8 


1.7 ±0.2 


Experiment, Dm 


[001] 


[111] 


Dqoi I Dill 


Ref.^ 


(x 10-12 cm2/s) 


7.1 ±0.5 


5.3 ±0.3 


1.34 ±0.12 


Ref4 (x 10-12 cmVs) 


6.4 ±0.9 


4.4 ±0.5 


1.45 ±0.26 


Experiment, D-^ 


[001] 


[111] 


DooiZ-Diii 


Rcf.* (x 10^12 cmVs) 


29 ±3 


33 ±4 


0.88 ±0.14 



u,q,l^k{u^q^l^k) 



+buqbkq { — 'iBuqlBkql + QBulqBqlk + SBuklBkql 
— 'iBkulBukl) + buqbkl {4:BkulBqku + ^BqakBukl 



'2Bi.„,iB, 



kqu^qkl 



'iBnk„.B 



qku^kql ) 



VI. NUMERICAL RESULTS 

Because Eqs. (EUl, (EH), and can be evalu- 

ated numerically only for finite lattice sizes, we use finite 
size scaling to extract the infinite lattice limit. The ap- 
proach to the infinite lattice value is expected to follow 
a power law. For example, if we approximate the sums 
by integrals in Eq. (j^ . 



D 



(0) 
M 



d:^rb{vYz^F{r). 



(41) 



<r<L 



As r cx), F — > const. Therefore, 



D 



(0) 
M 



const X 



const X 



r dr \ — \ r~ 



1 1 

a L 



(42) 



A least squares fit to a power law describes very well the 
scaling with finite lattice size. For diffusion of magnetiza- 
tion, we were able to vary the lattice size, in increments 
of 2 lattice sites, between 1 and 81 lattice sites on an 
edge. For energy diffusion, we studied lattices with up to 
27 sites on an edge. 

The constant F in Eqs. and H40|l was taken as 
the mean value of the Fij over three layers of nearest 
neighbors. This averaging procedure works well since 



contributions from far-apart indices are suppressed by 
the b factors in the summations. In dimensionless units, 
this gave 



(40) 



X F = 



0.48 ± 0.05, (001) direction, 
1.17 ± 0.14, (111) direction. 



(43) 



The error, taken to be one standard deviation from the 
mean, is 10% for the (001) direction and 12% for the (111) 
direction. Therefore, the error in F'^ is about three times 
as much, or between 30% and 40%. As the magnitude of 
the second order terms is between 10% and 20% that of 
the zero order terms, the overall error in approximating 
F„ « F is less than 10%. 

The results, including errors due to fitting and approx- 
imations, are summarized in Table Our magnetization 
diffusion coefficient agrees well with experimental values. 
However, the energy diffusion coefficient that we obtain 
is significantly smaller than observed experimentally, al- 
beit larger than the magnetization diffusion coefficient. 
The relative orientation dependence is within the experi- 
mental range for magnetization, but disagrees drastically 
for energy. 



VII. CONCLUSION 

We have presented a hydrodynamic approach to study 
the long- wavelength spin dynamics in a lattice of dipolar- 
coupled spins-1/2 in high magnetic field and at high tem- 
perature, motivated by recent experiments on coherent 
nuclear spin transport in calcium fluoride. The Kubo 
formula, Eq. (QJ, for the diffusion coefficients applies to 
the physical regime probed experimentally and rests on 
the assumption that the time evolution of the system 
is ergodic. We developed a perturbation theory for Eq. 
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that is equivalent to the approach of LGK but sim- 
pUfies the calculations enormously. This allowed us to 
obtain the diffusion coefficients for magnetization as well 
as energy to leading order in the flip-flop term of the 
Hamiltonian, and estimate the first pcrturbative correc- 
tion. The result for magnetization diffusion agrees with 
experiment to within its degree of accuracy. The result 
for inter-spin energy diffusion is larger than that for mag- 
netization diffusion, in qualitative agreement with the ex- 
periment, ft does not, however, describe the experiment 
quantitatively. 

One possible reason for the disparity is the ergodicity 
assumption, which is implicit in our choice of how to take 
the average in the correlation functions appearing in Eq. 
(0 . Our assumption of Eq. jSJ leads to the equivalence 
of the non-equilibrium average of the conserved density 
to its equilibrium correlation function, as is usual in lin- 
ear response theoryiii For magnetization diffusion, the 
low polarization results in a very low density of polarized 
spins in a completely randomly polarized background. 
The sparsity of polarized spins means that their effect on 
each other is negligible, and they can be treated as in- 
dependent. This suggests that statistical averaging over 
a complete infinite temperature ensemble is the correct 
procedure. For energy diffusion, the correlations inherent 
to the initial states used in the experiments of Boutis et 
al. may require the use of an ensemble that is a subset of 
the full Hilbert space. This would imply a modification 
of Eq.® for the density matrix at finite time, and needs 
to be investigated further. 

The most evident source of inaccuracy in our calcu- 
lation is our truncation of the perturbation series. Our 
estimate of the next-to-leading order correction to en- 
ergy diffusion does not rule out the possibility that a re- 
summation of our perturbation series would explain the 
experiment quantitatively. However, currently available 
non-perturbative (e.g. Bennet and Martini) and resum- 
mation (e.g. Borckmans and WalgraeSii) methods are 
much too cumbersome to treat the diffusion of inter-spin 
energy. They are also plagued by their own uncontrolled 
approximations, such as the ad-hoc replacement of cer- 
tain correlation functions by gaussians in order to sim- 
plify the calculations. This state of affairs, coupled with 
the great success of LGK at predicting the magnetiza- 
tion diffusion coefficient based on an expansion in the 
flip-flop term, and earlier that of Lowe and NorbergS& 
at fitting the shape of the FID by a similar expansion, 
motivated us to systematize their approach and apply it 
to inter-spin energy diffusion. For the final word on this 
problem, we await either a new non-perturbative method 
or a tractable technique for summing our perturbation 
series to all orders. 
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APPENDIX: EQUIVALENCE TO LGK METHOD 

LGK compute the time evolution of the average den- 
sity, (5(k, i)) = tr {^p(— k, t)5(k)}, starting from one of 
the non-equilibrium states, Eq. Q or Eq. Q). In general, 
we write 

Sp{t = 0) = e J d=^rcos(k • r)S'(r) = ^^(k) -I- ^S{-k). 

(A.l) 

The average of any Schrodinger operator A at time t is 



{Ait))= tr{5pit)A}. 
If {S{r,t)) satisfies a diffusion equation, 
d 



dt 



=I?V^(5(r,t)), 



(A.2) 



(A.3) 



then the time dependence of (5'(k, t)) [the Fourier trans- 
form of {S[v, t)) at wavevector k] is given by 



d_ 
dt 



(5(k,t)) 
{S{k,t)) 



-Dk^{S{k, t)) 
-'='^*(5(k,0)). 



(A.4) 



Taking the time derivative of the last equation and rear- 
ranging terms, we obtain LGK's expression for the diffu- 
sion coefficient, 



D = hm 



1 \ {S{k,t)) 



(A.5) 



The right hand side is in fact independent of time for 
t greater than the short timescale defined by the inverse 
spin-spin coupling. We keep the explicit time dependence 
to remind us that t is long, but finite, while k tends to 
zero. In other words, we take the t oo limit after the 
k —f limit. 

In the Heisenberg representation, we have 

{S{k,t)) = tr{e-'«*(5p(0)e'^*5(k)} 
= ^tr{S{k,-t)S{KO)} 

+ |tr{5(-k,-t)5(k,0)}. (A.6) 

Because of translational invariance, only the second term 
on the right side of Eq. ljA.6|) contributes, and we have 



(5(k,t)) = -tr{5(-k,-i)^(k,0)}. 



(A.7) 



We take the time derivative of Eq. (|A.7|I . and use 

dS{t)/dt^i[H,S{t)], to find 



(S(k,t)) 



-I- tr {e-*^* [n, Si-k, 0)]e''^*S{k, 0)} 
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■tr{S'(-k,0)5'(k, t)} 



Substituting Eqs. HA.9|) and (jA.10|) into Eq. HA.5|) gives 



dt' tr{S'(-k,0)5(k,i')} 



■tr{S'(-k,0)S'(k,0)}. 



(A., 



The second term on the right-hand side of the last 
equation vanishes. Substituting the continuity equation, 
^(k, t) +i'k- j(k, t) — 0, and taking k = kz, we obtain 



Ipdi' tra(k,^Oj,(-k,0)} 
-™ tr{5(k,0)5(-k,0)} ■ 



fc^O 



(A.11) 



fc2 / dt'tra(k,t')j.(-k,0)}. (A.9) 



(A.IO) 



(5(k,t)) 



According to Eq. I|A.6|) . we further have 



(5(k,0)) = -tr{5(k,0)5(-k,0)}. 



Taking the hmit i — > oo, we obtain the standard form of 
the Kubo formula, Eq. , which proves the equivalence 
of the two approaches. 
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